WHAT IS CLAIMED IS: 

1. A numerical calculation method for a physical quantity U practiced on a 
computer by solving A-U = f, wherein A is a coefficient matrix (in N rows by N columns; 
wherein N is a positive integer) obtained through discrete of a partial differential equation 
5 to be satisfied by said physical quantity U, and f is an inhomogeneous term (source term), 
comprising the processes of: 

setting an initial value of said physical quantity U; 

setting 0 as an initial value of a number m of repeating times, giving 0 as an initial 
value of a perturbation quantity ^and setting (f - A • U^) as an initial value r° of a 
10 residual r; and 

repeatedly executing a first step and a second step while incrementing said 
number m of repeating times until an approximate solution U*" is converged, 
wherein said first step includes the steps of: 

obtaining a predicted approximate value \|/" of A- = r"* through repeated 
15 calculation performed by a first calculation unit including an internal solver, and 
said second step includes the steps of: 

obtaining, from said predicted approximate value v|/", a corrected approximate 
value for minimizing norm of a residual r^ through an optimization routine 

performed by a second calculation unit; and 
20 giving (U*" + ^™) as an approximate solution U'"^^ and giving (r^ - A- ^J"") as a 

residual 

wherein in said second step, obtained elements of a vector sequence A-^*" are 
sampled by a given sampling method to be stored in a memory, and 

a residual minimization coefficient ai*" (wherein 1 = 1, L) used for obtaining 
25 said corrected approximate value ^"^ is approximately obtained by using elements of a 
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vector sequence A- (wherein k = m- L+ l, ...,m-l) stored in said memory. 

2. The numerical calculation method of Claim 1, 

wherein in sampling of elements bi, b2, ... and bw of said vector sequence A-^"" 
performed in said second step, elements bi (wherein ie fl) are selected, whereas a subset Q 
is defined as follows: 

Q={i:mod[i,lg] = l} U {i : |fi/aii| > p} 

wherein Ig is an integer, p is a real number, fj is an element of said source term 
and aii is a diagonal term on the ith row in the ith column of said matrix A. 

3. A numerical calculator for a physical quantity U by solving A-U = f, wherein 
A is a coefficient matrix (in N rows by N columns; wherein N is a positive integer) 
obtained through discrete of a partial differential equation to be satisfied by said physical 
quantity U, and f is an inhomogeneous term (source term), performing the processes of: 

setting an initial value U° of said physical quantity U; 

setting 0 as an initial value of a number m of repeating times, giving 0 as an initial 
value of a perturbation quantity ^and setting (f - A • U^) as an initial value r° of a 
residual r; and 

repeatedly executing a first step and a second step while incrementing said 
number m of repeating times until an approximate solution U*" is converged, 
wherein said first step includes the steps of: 

obtaining a predicted approximate value of A-^ = r"™ through repeated 
calculation performed by a first calculation unit including an internal solver, and 
said second step includes the steps of: 

obtaining, from said predicted approximate value v|/^, a corrected approximate 
value ^"^ for minimizing 1? norm of a residual r"™ through an optimization routine 
performed by a second calculation unit; and 
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giving (U"" + ^i"') as an approximate solution U'"^^ and giving (r"" - A-^*'') as a 
residual t'^''\ 

wherein in said second step, obtained elements of a vector sequence A-^"™ are 
sampled by a given sampling method to be stored in a memory, and 

a residual minimization coefficient ai"" (wherein 1 = 1, L) used for obtaining 
said corrected approximate value ^"^ is approximately obtained by using elements of a 
vector sequence A- (wherein k = m- L+l, ...,m-l) stored in said memory. 

4. The numerical calculator of Claim 3, 

wherein in sampling of elements bi, ba, ... and bn of said- vector sequence A-^"* 
performed in said second step, elements bi (wherein i e £2) are selected, whereas a subset Q 
is defined as follows: 

n={i:mod[i, lg]=l} U {i : |fi/aii| > p} 

wherein Ig is an integer, p is a real number, is an element of said source term 
and aii is a diagonal term on the ith row in the ith column of said matrix A. 

5. A recording medium that stores a numerical calculation program for a 
physical quantity U by allowing a computer to solve A-U = f, wherein A is a coefficient 
matrix (in N rows by N columns; wherein N is a positive integer) obtained through discrete 
of a partial differential equation to be satisfied by said physical quantity U, and f is an 
inhomogeneous term (source term), 

wherein said numerical calculation program makes said computer to execute the 
processes of: 

setting an initial value if of said physical quantity U; 

setting 0 as an initial value of a number m of repeating times, giving 0 as an initial 
value of a perturbation quantity ^ and setting (f - A • U°) as an initial value r° of a 
residual r; and 
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repeatedly executing a first step and a second step while incrementing said 
number m of repeating times until an approximate solution U"* is converged, 
wherein said first step includes the steps of: 

obtaining a predicted approximate value \|/" of A- ^ = r^ through repeated 
5 calculation performed by a first calculation unit including an internal solver, and 
said second step includes the steps of: 

obtaining, from said predicted approximate value v/", a corrected approximate 
value for minimizing I? norm of a residual r^ through an optimization routine 
performed by a second calculation unit; and 
10 giving (U™ + ^J™) as an approximate solution U*""^^ and giving (r^ - A- ^'") as a 

residual 

wherein in said second step, obtained elements of a vector sequence A-^*" are 
sampled by a given sampling method to be stored in a memory, and 

a residual minimization coefficient ai"^ (wherein 1 = 1, L) used for obtaining 
15 said corrected approximate value is approximately obtained by using elements of a 
vector sequence A- ^ (wherein k = m- L+l, ...,m-l) stored in said memory. 
6. The recording medium of Claim 5, 

wherein in sampling of elements bi, b2, ... and bw of said vector sequence A-^*" 
performed in said second step, elements bi (wherein i e £2) are selected, whereas a subset Q. 
20 is defined as follows: 

Q={i:mod[i, lg] = l} U {i : Ifi/anl > p} 

wherein Ig is an integer, P is a real number, fj is an element of said source term 
and aii is a diagonal term on the ith row in the ith column of said matrix A. 
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